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Abstract. In this paper, the author proposes a numerical method to solve a parabolic system of 
two quasilinear equations of nonlinear heat conduction with sources. The solution of this system 
may blow up in finite time. It is proved that the numerical solution also may blow up in finite time 
and an estimate of this time is obtained. The convergence of the scheme is obtained for particular 
values of the parameters. 



I. Introduction 

The purpose of this paper is to study the numerical behavior of the solution of a nonlinear 
reaction diffusion system with nonlinear source terms. 



Let Q a smooth bounded domain in IR d . We consider the system: 

u u - A< +1 = av? +1 , i6fi, t > 
v lt - A< +1 = a< +1 , x Gil, t>0 
ui(x, 0) = u w (x) > 0, x E Q 
t>i(x, 0) = v w (x) > 0, iffi 
Mi = v i = 0, x G d£l, t > 

with z/, /i > 0, a > 0. 

Samarskii and al. [5] have studied this system and obtained the following results: 

If Ai denotes the first eigenvalue of the Dirichlet problem: — Ap = Xp, x G f2, p = 0, x G dfl and if 

a > Ai, the problem has no global solutions and there exists To > such that 

, (\\u» + \t,.)\\* + H + \t,.)\\") = +oo 



In [3], [1], we have proposed a numerical method to solve a quasilinear parabolic equation with 
blow-up of the solution. The numerical solution is computed by using the function u = ; so 
the nonlinearity is reported on the derivative in time. This solution has the same properties as the 
exact solution, in particular blow-up in finite time. We generalize this method to the system of two 
equations. 



For what follows, it is more convenient to work with a transformed equation. Let u = u\ 



_ _ i 



in - ^pj, p = — j-j-; then we get m,p G]0, 1[ and we suppose that p < m (or p > v) 



l 
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Then (u, v) satisfies the following system: 



mii m_1 M[ + Au = av 
pv p ~ 1 Vt + Av = au 
u(x, 0) = u^fa) = u (x), 
v(x, 0) = Viq 1 (x) = v (x), igU 



where A is the operator —A of domain D(A) = Hq(Q) fl H 2 (Q). 

An outline of the paper is as follows: In Section 2, we study the asymptotic behavior of the solution. 
In Section 3, we define a numerical scheme and prove the existence of the solution of this scheme. 
The section 4 is devoted to the properties of the scheme, in particular, the existence of a numerical 
blow-up time in the case a > Ai- Finally, in Section 5, we study the particular case p = m and 
prove the convergence of the scheme in that case for a class of initial conditions. 



2. ASYMTOTIC BEHAVIOR OF THE SOLUTION 

Given u ,v G L°°(fl), a couple (w, v) is a weak solution of (11.11) on [0,T] if u, v G L°°((0, T) x Q) 
and 

/ (u m (f)t — uA(f) + av(f))dxdt = / u(x,t)(f)(x)dx — / uo(x)(p(x)dx 
o Jn Jn Jn 

/ (v p (f) t — vAcj) + aucj))dxdt = / v(x,t)(J)(x)dx — / v (x)(f>(x)dx 
'o in 

for all G L7 2 ((0,T) x Q) n C^QO, T] x = for x G <9R 

This problem admits a local solution and from the maximum principle, we get u(t),v(t) > in 
Q, < t < T. 

We prove the following results: 

l) if a > Ai, the solution blows up in finite time T; we get: 



l im [ ( JH!_ u m+i(t) + -JL^ V P+I( t )) dx = +00 



ll) if a < Ai, the problem has a global solution which tends to when t — >oo 

ill) if a = Ai, the problem has a global solution (u,v) which tends to Qp\ when t — >oo where p\ is 
the first eigenfunction of A (Api = \\pi and HpiH^^) = 1) and 6 is a constant depending on the 
initial condition. 



We introduce the functions <3>, Z defined on [0, T] by 
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$(t) = f ( -^— u m+ \t) + -^—v p+1 (t) J rfa; 



m— 1 



(2.1) = = ( / ( — ^u m+1 m + -^-w p+1 (t) ) x 

and the functional defined on Hq(Q) x ^(^by 



(2.2) J(u,v) = I (\Vu\ 2 + \Vv\ 2 -2auv)dx 

Jn 



Lemma 2.1. The function $ is convex and the function Z is concave 

Proof: We prove that the second derivative of $ is nonnegative. 
We have: $'(£) = f n (mu m u t + pv p v t )dx 

By multiplying the first equation of (11.11) by u, the second by v and integrating over Q, we 

$'(t) = / (mu m u t + pv p v t )dx = — I ((Au — av)u + (Av — au)v)dx = 
Jn Jn 



Then we deduce: 



/ (| Vu\ 2 + | Vf | 2 — 2auv)dx = — J(u, v ). 
Jn 



<&"(*) = -2 / ((Am - + (Aw - au)v t )dx 



n 



= 2 [ (—(Au-av) 2 u 1 - m + -(Av-au) 2 v 1 - p ) dx 
Jn\m p J 

and > 0, Vt > 0. The function $ is convex. 

Besides, we have: 

Z'(t) = I^l^i$m- 2 /(™+i)$'( t ) 
m + 1 

and Z"(f) = (^(t))- 2 /^ 1 )- 1 (^($'(t)) 2 - $(t)$"(t)) . 
By using Cauchy-Schwarz inequality, we get: 

(Au - av)udx < ( m + 1 / (Am - cro) V~ m cfa^) f / —^—u m+1 dx) 



m Jo J \ Jo m + 1 
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and ^ 

J (Av-au)vdx< {^-y 1 j (Av-aufv l ~ p dx^j -E— V p+l dx^j 

We deduce: 

($'(t)) 2 < ( [ (Au - avfu l - m dx + f (Av - aufv l ~ p dx\ $(t) 

V m Jn V Jn J 

and 

-2/(m+l) 



Z"(t) < 2-^^ ( I ( ^^u m+l + -^v p+1 )dx) I (Av- aufv^dx. 

~ (m+l) 2 p \J n \m+l p+1 ' ' 1 



v. 



Since p < m, the function Z" is nonpositive , that is the function Z is concave. 

In the two next propositions, we prove that in the case a > Xi, the solution (u, v) of (11. ip blows up 
in finite time and obtain estimates of this time for a class of initial conditions. 

Proposition 2.2. If a > X% and if the initial condition satisfies J(uq,vq) < 0, the solution blows 
up in finite time T such that 

m m+l i V p+1 l , 

Un H -f n dx 



(23) T< 1+m in V™+! ° P+l ° 



1 — m 



(|Vm | 2 + |Vt> | 2 - 2au v ) dx 



Proof: The function Z' is nonincreasing, so we get: 

(2.4) Z(t) - Z(0) < tZ'(0) 

-2 

1 — m 



and Z'{0) = — ($(0))m + l J(u ,v ). 

1 + m 

The inequality (12.41) may be written as 



n 



m— 1 m — 1 

m+l III 1YI T) \ \ 



m + l p+l / / \Jn \ m + 1 p+l 



+ _Ji_^ )P +i j dx j < ( / ( < +1 + -^v p Q +1 ) dx 
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If a > Ai , the set S = {(u ,v ) 6 Hq(£1)/ J(u ,v ) < 0} is not empty. If J(u ,v ) < 0, the 
right member of (12. 5p becomes equal to zero for a finite time and so the solution blows up in a finite 
time T such that 



T < 



m+1 
1 — m 



rn 



P 



,p+i 



m 



p+ 1 



dx 



-J(uo,v ) 



Proposition 2.3. : If a > \±, the solution blows up infinite time T and satisfies the inequality: 



1 l-m 1 

Proof: If J(wo, fo) < 0, from Proposition 2.3, we know that the solution blows up in finite time. If 
J(uq, vq) > 0, the result is obtained by using the same proof as Friedman-McLeod in pp. 

Since the function Z is concave, it satisfies 

Z(t) > ^^Z(0) + -Z(s), < t < s < T and lim Z(s) = 0. 
s s s — >T_ 

We deduce 

m— 1 m— 1 

that is the inequality (12. 3p . 

Proposition 2.4. If a < X\, the solution tends to u>/ien £ — > + oo. 

Proof: If Ai > a, we get: J(u, v) > (Ai — a) J Q (u 2 + v 2 )dx > and < 0. Since the function 
$ is convex, the function $' is nondecreasing and lim <&'{t) = I < 0. 

t >+oo 

We deduce: $(t) < $(0) + /£, Vt > 0; since > 0, we obtain that I = and 
(2.6) lim J(u(t),v(t)) = 0. 

t — >+oo 

If Ax — a > 0, then, we get ^ lim J n (w 2 + v 2 )dx = and the solution tends to when t — > + oo. 
Proposition 2.5. If \\ = a, then lim u(t) = lim v(t) = Qp\ in L 2 {Vt) where 9 is a constant 

t >+oo t >+oo 

depending on the initial conditions. 

Proof: If Ai = a, we get from (12.61) that lim J(u(t),v(t)) = and is bounded for t > 0. 

t >+oo 

By interpolation, we obtain 



6 
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/ u 2 dx<5 [ \Vu\ 2 dx + d(6) ( f u m+1 dx 
Jn Jn \Jn 



2/(m+l) 



r r / r \ 2/(p+1) 

/ v 2 dx<5 / |Vv| 2 cfe + C 2 (5) / v p+1 cfe 

where Ci(5) and C 2 (S) are constants depending on fl and m and p respectively and we get: 

r / r \2/(m+l) , . \ 2/(p+l) 

(u 2 + v 2 )c£r < <JJ(«(t), v(t))+2<JAj / uvdx +d(S) I / u m+1 cix +C 2 (5) / v p+1 dx) 

n Jn \Jn J \Ju J 

We deduce 

, , , x 2/(m+l) x . \ 2/0+1) 

(1 - <5Ai) / (m 2 + v 2 )cfe < 5J(u(t), v(t)) + d(5) / u m+1 dx + C 2 (S) / v p+1 dx 



Jn \Jn / vn 

Since J(u(t),v(t)) is bounded, if 5 is chosen such that 1 — 8Xi > 0, we deduce that u{t) and v(t) 
are uniformly bounded in L 2 (fl) and in #0(0). So, we can extract subsequences t n — > + oo such 
that u(t n ) and v(t n ) converge weakly in Hq(Q) and strongly in L 2 (fl) to Z\ and z 2 respectively. 

We have : J(z 1 ,z 2 ) < lim J(u(t), v(t)) = 0, that is j n (|V.zi| 2 + | V^ 2 | 2 — 2\ 1 z 1 z 2 ) dx = or 

/ (\V zi\ 2 — Xizl) dx + / (| Vz 2 \ 2 — X\zl) dx + Ai / (z\ — z 2 ) 2 dx = 

Jn Jn 

We deduce that z\ = z 2 = Qp x - 

By multiplying the two equations of (11.11) by p 1 and integrating over fl, we get: 
dt 



J {u m + v^pxdx^j = and then 



u m (t) +v p (t)) Pl dx 
Hence we obtain 



n 



/ ( u o + v^) Pl dx, for all t > 0. 
Jn 



e m p rn+l + QP pP+^ dx = / + ^^fa 

n Jn 
and there exists a unique positive value of 6 satisfying this equation; we deduce the proposition 



3. Definition of a numerical scheme 

The classical Euler scheme cannot blow up in a finite time, so we generalize here to a system the 

numerical scheme used in [3]. 

The first equation of (11.11) may be written: 

m _ u m- P+ id , u = av 



1 — p dt 
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- — v—(v p l ) + Av = au. 



1 — p dt 

So, we discretize the two derivatives in time in the same manner: 

If (u n ,v n ) is the approximate solution at the time level t n = nAt, ( where At is the time step), the 
approximate solution at the time level t n+ i is solution of the system: 

(3.1) -^u n + x u™- v {v?- x - u p - + \) + AtAu n+1 = aAtv n+1 , 

1 — p 



(3.2) i~p Vn+1 ^ 1 1 ~~ + ^ tAv n+i = cnAtu n+1 . 

We prove first that the system (13.11) . (13. 2p has a unique positive solution if u n ,v n are positive in Q. 
We need several lemmas. For what follows, we denote: ||t>|| r = IHI^n) • 

Lemma 3.1. If the functions u n andv n are positive on Q, continuous in Q and satisfy the condition: 

In o\ II 111— m n l|l-P „ m P 

(3-3) \\Un\\oo IHloo < 



a 2 (l -p) 2 At 2 

then the system $3. $3. ~E) has a positive solution u n+ i,v n+ i G C 2 (fl) 
Proof: Consider the functional defined on Hq(Q) x Hq(Q) by: 

(3.4) Jn(u,v)= [ (|Vd 2 + \Vv\ 2 -2auv)dx + 1 —— [ (mu^u 2 + pvET l v 2 )dx 

Jn {±-P)AtJ n 

and let us denote K = j (u,v) G (H^Q)) 2 / J {mu™~ v \u\ p+1 + p\v\ p+l )dx = lj . 

Since (u n ,v n ) satisfies (13. 3L we get : mii^S 2 + pv p ~ 1 v 2 — 2a(l —p)Atuv > and 
J n {u,v)>0, J n {\u\ , \v\) < J n (u, v), Vu, V G Hq(Q). 
Now, we consider the problem: 

(3.5) min J n (u, v) 

From the preceding remark, the solution of (13.51) . if it exists, will be positive. 

We denote <p(u, v) = J n (mu™~ p u p+1 + pv p+1 )dx. 

If (u,v) is a solution of problem (I3.5p . there exists A G R such that: 

' dJ n (u,v) x 50 
<f>(u,v) = 1 
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Then, we get for any ip e Hq(Q), 



f TTl f 



/ (Vt> Vip — auip)dx + — 
Jn (1 



V 



p)At 
and 

(mu p ~ l u p+l + pv p+1 )dx = 1. 



u™~ v uipdx = Xm(p + 1) / u™~ p u p i>dx : 
n if7 



f,^ l vil)dx = Xp(p + 1) / v p ipdx 

n 



Hence, (w, -0) satisfies the equalities: 

777 

Au-av+ = Am(p + 1)<^ P , 

AO — mi + — — f = Xp(p + 

(1 -p)At 

and (m, -0) G C 2 (fi). Besides J n (u,v) = X(p + 1). 

A solution of (13.11) . (I3.2p is then defined by u n+ \ = ju, v n+ \ = 7t>; we get 

Au n+l - av n+l + ——— u ™- l u n+1 = X^-Pmip + l)CX+i. 



l-p)At V " lvn+1 = At * Pp<yP + 



Au n+ i - au n+ i + — ZT7n v n v n+i = A7 p p(p + l)u 

and 

(3.6) 7 = 



x \ l/(l-p) 



(1 - p)AtJ n (u,v) 
Hence, the numerical scheme admits at least one positive solution. 

Before proving the uniqueness of the solution, we first prove the existence of bounded supersolutions 
and of maximal solutions. 

Lemma 3.2. // the hypotheses of the lemma (E2J) are satisfied, the system ( Iff. II) . (EHj) admits a 
constant super solution. 

Proof: Let {C^Cl) e IR+; (C*, C 2 ) will be a supersolution of the system (IQl if these 

constants satisfy the inequalities: 

777/ 777 
(3-7) - ^ _ p)At ( C nY<~ P + ^_ p)M C n U n~ l > <*(%, 
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The first inequality may be written: 



We note x — — 



The second inequality may be written: 



[l-p)At n ~ (1 -p)At 

m 



hence this inequality may be satisfied only for x < At ||w n ||ro 1 



V „_i ol ^ p 



r 1 - - > TT-^TTT^r 1 



(1 -p)At n x - (l-p)At 

and may be satisfied only if — < — — ||u n ||*'7 1 . 

3 3 x (I -p)aAt 11 "°° 

So a necessary condition to obtain inequalities (13. 7p . (I3.8P is that 
(1 -p)aAt |L _ Ml _p m ! 



t; «lloo P < 71 t — t— , that is the condition (13.31) 



p 11 " llo ° (l-p)aAt 

(l-p)aAt nl-p. . m I, iim-1 

and \\v n I < x < — — \\u n I 

p 11 ra|lo ° (l-p)aAr 1 "°° 



If we choose = — — 7771 — r, then we get : 

" 'l ^-P) a, I, |,l-mV /(1 - p) 

1 Arx||u n 

m 

q2 _ X ll M "llno 



It remains to prove that x may be chosen in the interval [ ||w 



such that Cl satisfies : (C^f- 1 < v v ~ x - -- — ^Ai. 

x p 

In order to obtain this inequality, the parameter x must satisfy: 



(3-9) rT77T^ < "''"""^ 



1 — p)aAt m 
p l|Wn|lo ° ' (l-p)aAt 



V(l-P) - A ..I-pa^-H |il-m\V(l-p) 

1 1 ) 



Let us denote a = — -Ata \\v n \\~ p , b = 

p II "lloo ' (1— p)o 



I Nm-l, 



we have : a < x < b and the condition (13.91) may be written: 



V 



6 (l-p)/(l-m) m ] (l - 

l-p)aAt V 6/ ~ V(! -p)aAtJ \ xJ 



(i_p)/(i_ m ) 
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If we define the function / by 

/ \ (l-p)/(l-m) 

Ux) = P-—ab (m - p)/[1 - m) x p (x — b) + I - m . (x - a) 

JK 1 (l-p)aAt v ; \(l-p)aAtJ y 1 

the condition (13. 9p becomes: f(x) > 0. 

The function / satisfies f(a) < and f{b) > 0; so there exists x &]a,b[ such that f(x ) = 
and the couple 

/ \ 

I 1 1 ~s 1 1 V n 1 1 



^1 _ II "nlloo ^y; 

" ~~ / ._>, \ i/u-pV 5 



is a supersolution of the system (13. ip . (13. 2p . 



Lemma 3.3. System $3. II) . ([Of /ias a maximal solution (u, v) and any solution (u, v) satisfies: 
< u <u, < v <v. 

Proof: We use the same method as Keller in [2]. We consider the sequences defined by: n n+10 = 

Vn+1,0 = Cni 

. yy^ ^ 771/ p yjy ^ 

Au n +\j+\ + ^ _ Un+1 >i +1 = aVn + 1 J + (1 — p)At Un+1 ' jUn ' 



We get: 

rn rn 

A(u n+ i,i - u n+1 , ) + (1 _ p)At <~ 1 K+i,i - Wn+i,o) = «^ + _ p)At C" P {(ClY ~ < C„) ■ 

The second member of this equality is negative; we deduce from the maximum principle that: 
Wn+i,i < Mn+i,o- in the same manner, we get: f n +i,i < ^n+i,o- We prove recurently that the 
sequences (^n+ij)j>o an d ( ,y n+i,j))i>o are decreasing; in fact, we have: 

Tfl y 

aK+l,j - Un+lJ-l) + ^ _ p \^ t U n~ P ( U n+l,j ~ 

and the second member is negative from the recurrence hypothesis. 
We deduce: Wn+ij+i < Un+ij- Similarly, we get: f n +i,j+i < v n+l j. 

Since the two sequences (u n +i,j)j>o an d (v n +i,j)j>o are nonnegative, they converge to u and v and 
taking the limit when j — > + oo, we obtain: 



NUMERICAL SOLUTION OF A PARABOLIC SYSTEM WITH BLOW-UP OF THE SOLUTION 



11 



P J) 

(1 -p)At n (1 -p)At 

So, the functions u and v are solutions of the system (13. lft . (13.21) . 

It remains to prove that any solution (u, v) satisfies: < u <u, < v <v. 

Let (u,v) a solution of system (13.11) . (13 ,2p . we have: < u < C*, < v < C\\ we have the 
equalities: 

TTT TTl 
A(U - U n+l!j+ l) + _ p -} At U n~ 1 ( U - Wn+lj+l) = a ( V - V n+l,j) + (1 - p) At ^^ ^ ~ U ™ +1 ^ 

P P 

A{v - v n+1)j+1 ) + —— vP-\ v - v n+1)j+1 ) = a(u - u n+lij ) + ——— (vP - v p n+lj ). 

For j = 0, the second member of these inequalities is negative; then we get u < w n +i,i > v < u n +i,i 
and recurrently, we obtain u < u n+ xj, v < v n+ ij for any j > 0. It results: u <u, v < v. 

Theorem 3.4. If the functions u n and v n are positive, continuous in Q and satisfy the condi- 
tion ( Iff. 3\) , then system ( Iff. ( fff.ffj) has a unique positive solution. 

Proof: From the previous lemmas, we know that the system admits at least one positive solution 
and that any solution (u, v) satisfies < u <u, < v <v. 
We get: 

TTl f f 

(Au u — Au u)dx = = — — / u!^~ v u u{uP~ x — vP~ )dx + a I (vu — vu)dx. 

' (l-p)AtJ u 1 J Jn ' 

Similarly, we have: 



n 



P 



v v (v p 1 — If 1 )dx + a (uv — uv)dx = 
n Jn 



(1 -p)At 

We deduce from these equalities that J n (uv — vu)dx = and then u —u, v — v. 
Theorem 3.5. The numerical solution exists at least until the time 



I m \m-l P \P— 1 

Tt = mm — -A , — r\ p 

\a{l-p) a(l-p) 



with A = max (||w |L , IK'o 



Proof: We prove recurently that the solution (u n ,v n ) satisfy the inequality: Hindoo ? H^nlloo < 

a(l-p) xm _! «(l-p) X p-A x A 



»i/(i-p)" 



where n is defined by O — max 5 \ 

\ "' P J (1-Woj 

If this inequality is satisfied at the time level t n = nAt, if £ n +i0o ^ 1, the inequality ( 13. 3p is verified 
and the solution exists at the time level t n+ \. 
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The quantity 4>n+i wm be a supersolution of the system (13.11) . (13. 2p . if we have the two inequalities: 



m 



in 



[i- P )Ae n+1 n + (i- P )At 



This may be written: 



(f>l+i > max 



|l-p 



1 — Oi 



(1-p) 



m 



a i 1 1 ii 1— m 1 I x /"V a j- 1 1 ll 1— P 

At||w n |L I- a At k) n I 



But, from the recurrence hypothese, we get 

ii-P 



7> 



< 



m 



1 1— m 
'>i\\oo 



A, 



1— m 




m (i _t n Q )(i-™)/(i-p) 



and it is easy to see that this quantitiy is bounded by 4> n +i- 

II II 1—33 
\\ v n\\ 

In an analogous manner, we obtain that = — is bounded by <f) n +i. 

1 1 ~~ P A* II II 1-P 

1 - a At \\v n \\ y 

p 

So, the solution at the time level t n+ \ satisfies: H^n+illoo ? Il ?; n+i|loo — ^n+i and the numerical 
solution exists exists during a positive time interval. 



4. Properties of the numerical scheme 

In this section, we prove that if a > Ai, the numerical solution blows up in finite time. 
We define the functional ip n and F n by: 

V(p+i) 

if> n (u,v)= ( / {mu^- p u p+1 + pv p+1 )dx 



and 

J(u, v) 
[il) n {u,v)) 2 

Lemma 4.1. The sequence (F n (u n ,v n )) n > is nonincreasing. 

Proof: Since u n+1 = and v n+1 = yv, (u, v) G K, we get ip n (u n+ i, v n+1 ) = 7 and F n (u n+1 , v n+1 ) 
J(u, v). 

Besides from (13. 4p . we have J(u,v) = J n (u,v) — ^ — - ^ J (mu™~ 1 u 2 + pv^v^dx. 
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Hence, we get : J(u, v) < 



In 



^l(u n ,v n ) (l-p)At i/jl(u n+1 ,v n+1 ) 



(ran™ l u 2 n+l +pv p n 1 v 2 n+l )dx 



In addition, we have the equality: 

J n (u n ,V n ) = J(u n ,V n ) + — -^(ipniu^Vn))^ 1 



We deduce: 



(1 -p)At 



J{U,V) < —tt, r + 



(ip n (u n ,v n )) 



p-l Jn 



(mu™ l u 2 n+1 +pv p n l v 2 n+1 )dx 



\ 



\ 



ij 2 (u n+1 ,v n+1 ) 



By the Holder inequality, we have at once: 

\(P+l)/2 



uT P ul + + \dx < 



I <X\dx <([ v 2 n+1 v^dx 
Jn \Jn 



U n+l U n dx J 
! / 

(P+l)/2 



u™ +1 dx 



(1-p)/2 



Hence, we get: 



/ /• i \(P+l)/2 

and 



~ l dx 



vl^v p - l dx\ { I v p+1 dx 



v p+1 dx 



(P+l)/2 



(1-p)/2 



(1-p)/2 



< +1 <£r 



(1-p)/2 



(4.1) Vn(«n+l 5 «n+l) < ( m / u l+l U n 1 dx + p / X dx ) Vn "("n,^)- 

\ Jn Jn J 

We deduce: J(u, 0) < J ^ n,Vn> ) ; that is F n (u n+ i, t>„+i) < F n (u n ,t> n ). 

Lemma 4.2. For n > ; we /iave i/ie estimate: 

(4.2) (l - p)AtF n (u n+1 , v n+ i) < ^ _1 (m„ + i, f n +i) - C 1 ^, ^n) < (1 - p)AtF n (u n , v n ). 



Proof: i) We prove first the right inequality. We have: i/;(u n+ i, f„+i) = 7; from (13.61) . we obtain 

i/j n (u n+ i,v n+1 ) = ((1 - p)AtJ n (u, {;)) 1/(p ~ 1} 



and 



^ u n +i) <(l-p)At 



i/;P- 1 (u n+1 ,v n+1 ) < (l-p)At \^F n (u n ,v n ) + (1 4)At ^n~ 1 ( u »» t, »)) • 



that is 
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a)Multiplying (13.1 p by u n+ i and (13.21) by v n+ \ and integrating on fl, we get: 

/ u 2 nJrl u™~ l dx - m \ u p n X\u™~ p dx+ / (\Vu n+1 \ 2 - au n+1 v n+1 ) dx 
Jn I 1 ~ P)^t Jn Jn 



(l-p)At 



(1 -p)At 
Hence , we get : 



P J v 2 n+1 v p l dx - _^ At J v p n X\dx + (| Vv n+1 \ 2 - au n+1 v n+1 ) dx = 



,. , {mu 2 n +1 u™ 1 +pv 2 n +1 v p l )dx-- — l —— ^J p+1 {u n+1 ,v n+1 ) + J{u n+1 ,v n+1 ) 

(l-p)AtJ n (l-p)At 

By using (14. ip . we deduce: 

CVn^'n) - 1p P ~ 1 (u n+1 ,V n+1 ) + F(u n+1 ,V n+ l) < 0. 

This concludes the proof. 

Lemma 4.3. The sequence {J{u n ,v n )) n >o is nonincreasing 
Proof: In [3] , we have proved the inequality: 

Va, b G R + , a p -\b - a) 2 < a p+1 - b p+1 - I^fc 2 ^- 1 - a p ~^. 

1 — p 

We deduce from this inequality: 

/ {u n+1 - u n fuZ~ x dx < [ u^dx- [ u p n + + \u^ p dx-\±P [ ul^u^-ul-^uZ 
Jn Jn Jn i — p Jq 

and 

/ (v n+1 -v n ) 2 v p ~ l dx< [ v^dx- [ v£\dx-\±£ [ v^vl'l-v^dx. 
Jn Jn Jn i — p Jq 

Since u n+ i and v n+ i are solutions of (13.11) . (13.21) . we get: 

/ M n+i(«n+i- M r 1 )«r P ^ = - — / (|Vu n+ i| 2 - au n+1 v n+1 )dx 



and 



~ v n l )dx = — — [ (\Vv n+1 \ 2 - au n+l v n+l )dx. 
V Jn 
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So, we obtain: 

m / (u n+ i - u n ) 2 u™- l dx + p (v n+ i - v n ) 2 v p ~ x dx 
Jn Jn 

< ip p+1 (u n ,v n ) - ip p+1 (u n+1 ,v n+1 ) - (1 +p)AtJ(u n+ i,v n+ i). 

From the inequality [3J :Va, b G M + , a p+1 - b p+1 < ±±| a 2 ^ 1 - a?" 1 ), 
we deduce: 

ip p+1 (u n ,v n ) -tp p+1 (u n+1 ,v n+1 ) < -^-i)l(u n , w„)« _1 (m„ + i,w„ + i) - ij p -\u n ,v n )) 

1 — p 

and by using (14. 2L we get: 

(4.3) tP P+ \Un,V n ) - Vn\u n+1 ,V n+l ) < (1 + P )AtJ(u n , V n ) . 

So, we get 

m / (u n+1 - u n ) 2 u™~ l dx +p / (v n+1 - v n ) 2 v p n ~ l dx < (1 + p)At (J(u n ,v n ) - J(u n+1 ,v n+1 )) . 
Jn Jn 

We deduce: J(u n ,v n ) > J(w„ +1 ,t)„ +1 ). 



We note 
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Lemma 4.4. // J(u ,t>o) < 0, t/ie sequence ($„)„>o is increasing. 
Proof: We have the equality : 

(4.4) e 1 ^, «») = / (jn< +1 + K +1 )^ = (p + l)*n + / u^dx. 

Jn m + 1 J n 

Besides, we get: 

(4.5) V n +1 (u n+1 , v n+1 ) = [ (mu™-" p < + + \ + Pvl + + \)dx <(p+ l)$ n+1 + m(m ~ P) / uT'dx. 
and we deduce: 

$n - $n+l < —— (^n +1 K, «n) ~ ^n+l)) • 

P + 1 

By using ( 14. 3 j) , we obtain $ n — $, t +i < AtJ{u n ,v n ). 

If J(mo,' i; o) < 0, since the sequence (J(u n , f n ))n>o is nonincreasing, we deduce that the sequence 
($ n )„> is increasing. 



16 MARIE-NOELLE LEROUX 

Lemma 4.5. If J(u ,v ) < 0, for n > 0, we have the inequality: 

(4.6) (^/^ - ^)/C^D) < ^V^n, t*0 «~ W V n+1 ) - ^(U^ V n )) 

Proof: This inequality may be written: 

Vn + \u n ,v n ) + (m + p +1 )< p + - 1)/(p+1) < ^K,^)C~Wi,^+i) + (m + 1)$„. 
By using (14.4D and (14.51) . we obtain that a sufficient condition to satify this inequality is: 

(p + 1)$„ + (m - pK + (m + l)$2/(f +1 )$( p + - 1)/(p+1) 

< ((p + l)$ n + (m - pK) 2/(p+1) ((p + l)$ n+ i + (m - p)/i„) (p - 1)/(p+1) + (m + l)$ n 
171 1 ,m+l. 



with u n = / u™ + dx. 

m + 1 J n 

Tt Tf \ ^ n *-\. \ ■ ■ ■ + (p+ l)$ n + (m-p)/i n $ n 

It J(w ,Vo) < U, the sequence { < P n )n>o is increasing, so we get: — > . 

(p + lj$n+i + {m -p)n n $ n+ i 

Hence, it is sufficient to prove: 

(p + 1)$„ + (m - p)/i n + (m + IJ^VCp+iJ^D/Cp+D 

< ((p + l)$ n + (m - p)/i„)) $(f-i 1)/(P+1) *£ 1_P)/(1+?,) + (rn + 1)S„ 
that is ^(^T^ - $r p)/(1+p) ) < *n(^f )/il+p) ~ &t pmi+p) ) 
and this inequality is satisfied since \i n < $ n . 

Lemma 4.6. If J(uq,vq) < 0, we have the estimate: 

(4.7) $f p+1 > (sfe 1 ^ 1 ' - $t 1)/(p+1) ) < ^AtJ(u n , v n ) 

Proof: We deduce the estimate immediately from (14.21) and ( 14. 6ft 

Theorem 4.7. // J(ito, fo) < 0, the numerical solution blows up in a finite time such that 

1 + 171 Jn\m+l U + p+l V ) aX 

l-p -J(uq,v ) 
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Proof: From g^D, we get $^ )/(p+1) < ^- 1)/{p+1) + ^At$" 2/(p+1) J{u n , v n ) 

and since the sequence (J(u n , v n )) n > is decreasing and the sequence (<3> n ) increasing, we get : 

and we deduce the estimate. 



Remark 4.8. In the case p = m, we obtain the same bound for the numerical blow-up time and 
for the blow-up time of the exact solution. In the case p < m, the bound obtained for the numerical 
blow-up time is inferior to the estimate obtained for the blow-up time of the exact solution 



5. The case p = m 

In the case p = m, the functionals ip n and F n are independent of n. We shall note them respectively 
ip and F: 

i J(u v) 

ij,(u, v ) = ( L m{u m+1 + v m+1 )dx) m+1 and F(u, v) = v ' \ . 

ip (u,v) 

Besides, we get: ip m+l {u n ,v n ) — (m+ l)$ n , n > 0. 

The estimate of the numerical blow-up is in that case the same as the estimate of the exact blow-up. 
5.1. Properties of the scheme. 

Proposition 5.1. If a > \\ and T* is the numerical blow-up time, we get the estimate: 
Proof: The estimate (14.21) may be written in this case: 

(1 - m)AtF(u n+ i,v n+1 ) < ^ m_1 (M n+ i,t; n+ i) - ip" 1 ' 1 (u n , v n ) < (1 - m)AtF(u n ,v n ). 

We deduce, since (F(u n , v n )) n >o is a nonincreasing sequence: 

^- l {u n+h v n+ j) - il)™- l (u n ,v n ) < (1 - m)jAtF(u n ,v n ), Vj > 0, 

il)™- 1 (u n -. i ,v n - i ) - ^^(u^Vn) < -(1 -m)iAtF(u n ,v n ), Vi > 0. 
So, we obtain: 

If T* = NAt is the numerical blow-up time, by choosing i = n, j = N — n, we get : 
(iV-n)V> m ->Q,0 < N4) l ~ m (u n ,v n ) or, ^(u n ,v n ) < -^-^""K,^). 
This is the same estimate as for the exact solution. 
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Proposition 5.2. If a < Ax, the numerical solution tends to when t — > + oo. 
Proof: In that case, we get: J(u,v) > (Ai — at) J n (m 2 + f 2 )<ix > 0. 

2 

Besides, since m < 1, for any G £ 2 (fi), we have: J n 2 dx > C(fi) (J n m+1 Gb) m+1 

and by using Young inequality, we get 

2 

2 2 



1 - m , „ x m + 1 



n 



2 m + 1 ( / (n m+1 + t; m+1 )rfa; ) < ( / u m+1 dx) m + 1 + ( / W m+1 da^ m + 1 



So, we get : F(u,v) > ^ " 



Besides, from ( 14. 2p . we obtain 1 (w n+1 , u n+1 ) > ?/> m 1 (-u n ,t>„) > ^ ^ At and 

O (i ZJ 

We deduce: lim L (u m+l + v m+l )dx = 0. 

n >+oo Jil 

Proposition 5.3. If a = Ai, then lim w n = lim v n = 9p\ where 9 is depending on the initial 

n — >+oo n >+oo 

conditions. 

Proof: We have the equality: $ n+ i - $ n = ^ (^> m+1 (it n+ i, v n +i) - ^ m+1 (M n , u n )) 

1 _|_ rji 

and from the inequality Vo, b G M + , a m+1 - 6 m+1 < -^— a 2 (b 711 ' 1 - a" 1 ' 1 ), 

1 — m 

we obtain: ^ 

-$„, < ^ 2 (M n+ i,W„ + i) Up m ~ l (u n) V n ) - V m-1 (l*n+l,Un+l)) < ~ At J(w n+ i , V n+ i) . 

I — m 

In the same manner as in proposition (12.5ft . we deduce that lim w n = lim v n = 9p\. 

n >+oo n >+oo 

We now obtain an estimate for 9 : 

By multiplying the equations (13. ip . (13. 2p by p\ and integrating on Q, we get: 

m f m f f 

- / u n+1 u™~ p x dx - - / Un+ipidx + AiAt / {u n+1 - v n+1 )p x dx = 0, 

1 — m J n 1 — m J n J n 

in i vn / / 

v n+1 v™~ p x dx - - / v™ +1 pidx + XiAt / (v n+1 - u n+1 )dx = 0. 



By using the inequality: u™ +1 < (1 — m)u™ + m« n+1 u™ 1 and the same inequality applied to the 
function v, we get 



/ u™ +1 p x dx < [ u™p x dx + At— / (v n+1 - u n+1 )pidx, 
Jn Jn m Jn 

\ v™ +1 pidx < \ v™pidx + At— / (u n+1 - Vn+Jp^x, 
Jn Jn m Jn 
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and then: / (u™ +1 + v™ +1 )pidx < 

Jn Jn 



u™ + v^)pidx, n > 0. 



We deduce: 6 m < 



n 



2 / p™ +1 dx 
n 



5.2. Convergence of the scheme. 

In this section, we obtain estimates on the numerical solution so we can extract by compactness a 
convergent subsequence. In order to prove that the limit is solution of the system (11. II) . we need 
an hypothesis on the initial condition, (this is due to the fact that we have a negative power in 
the scheme and we may not use a Holder inequality). If the initial condition does not satisfy the 
hypothesis, we observe numerically that this hypothesis is satisfied after a few times steps and the 
scheme again converges. 

Let us denote by 7T = inf TJAt) if T* is the existence time of the numerical solution. It 

0<At<At o 

follows from theorem [33] that T* > T\ > 0. Let T G [0,T*[, (T = AT At). We denote u At and v At 
the approximation of u and v defined by: 



"a/(/) = k + ^«Vi-<: 



l/m 



1/m 



Theorem 5.4. The sequences (um) an d (v At) ore uniformly bounded in C(0, T; Hq(Q)) and -f^ 1 (0, T; L 2 



Proof: Since T < T*, the functions (u n ) n >o and (f n )n>o are uniformly bounded in C(0, T;Q) 
and since J(u n ,v n ) is nonincreasing, we get: J(u n ,v n ) < J(uq,vq). We deduce that the sequences 
(Vw n ) n >o and (Vi> n ) n >o are uniformly bounded in L 2 (Q). 

We prove now that the sequences (— f±) and (^f ) are uniformly bounded in L 2 (0, T; L 2 (Q)). 
We have: 



du 



At 



dt 



L 2 (0,T;L 2 (Q)) 



JV-l 

E 

n=0 



tn + l 



n m 



2 M At 



2(l-m) 



m _ m 

At 



dxdt 



N-l 



n+1 ^ + U n ^ ) ( M n+1 M n ) <^ X - 



n=0 
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From the following inequalities, u l n ^ \u™ +1 — u™\ < \u n+ \ — u n \ and u\ m \u™ +1 — u™\ < \u n+ i — u r 
we deduce: 



dUAt 



dt 



L 2 (0,T;L 2 (fi)) 



< 



m 



1 JV ~ 1 f 

F— V / (u n+ i - u n fdx 

At n=0 J V 



< 



*a7 



m 2 At 



n=0 



.m— 1 / 



[Un+1 



\u 



1 1— m 



"Moo 



and we obtain analogous inequality for da*. 
Hence, we get: 



dUAt 



dt 



+ 



L 2 (0,T;L 2 (Q) 



At 



L 2 (0,T;L 2 (n)) 



< C (J(uo,Uq) - J(u N ,v N )) 



since Mr 



J(Kn<7V 



and (\\v„ 



)0<n<N 



are uniformly bounded. 



If a < Ai, we have: J(un,vn) > 0, 

if a > Ai, we get J(un,Vn) > (Ai — a) L (w^ + v 2 N )dx and this quantity is bounded from below; 
we deduce that the sequences are uniformly bounded in L 2 (0,T; L 2 (f2)). 

Since the sequences («At) and (t>At) are uniformly bounded in C(0, T; iJg (ft)) PI H x (0, T; L 2 (f2)), 
we can extract subsequences which converge to functions w and v in C(0,T;f2) if d = 1 and in 
(7(0, T; L r (fi)), (r < 2d/(d - 2) if d > 2, r = oo if d = 2,( Simon ®). 

In order to prove that the limits w, v are solutions of the system (11.11) . we use the same proof as 
in [3]; so we need to estimate the quantities 
is the object of the two next lemmas. 



^X" 1 ! and Ifer-^K 1 " 1 !- This 



Lemma 5.5. For n > 0, we have the inequalities: 



(5.1) 



Proof: These inequalities are proved recurrently; it is true for n — 0. Assume it is true at the order 
n — 1, that is : 



l—m 



l/(l-m) / \l/(l-m) 



The functions ( t^ 2 - ) w n , ( 7^ ) f n will be subsolutions of (11.11) if : 

\ fn+1 / \ Cn + l / 

<* + AtAu n - aAtv n < 0, 



+1 

m At 



1 — m t n 



m At 



v™ + AtAv n - aAtu n < 0. 



1 — m t n n 

But u n , v n are the solutions at the time level t n , so we get: 
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m At u-+At(Au n -av n ) = ( ^< + u n u-l - uA = ( u^u^ - ^ 



I — m t n 1 — m \t n J 1 — m \ t 1 

and --^—^v™ + At(Av n - au n ) = --^-v^v^v^ 1 - 

1 — m t n 1 — m tn 

From the recurrence hypotheses, the second members of these two inequalities are non positive and 

the inequalities (15 .ip are satisfied at the order n. 

Lemma 5.6. Assume that the initial conditions uq,vq satisfy: 
There exists a constant Cq such that 

(5.2) Au - av + C w™ > 0, Av - au + C v™ > 

Then, we have the estimates: 

rp , \ l/(l-m) / rp , \ l/(l-m) 



(5-3) u n+ i < u n , v n+ i < 

\ 1 2 — tn+1 / V 1 2 — l n+l 

with T 2 



(l-m)Co 



Proof: This lemma is proved recursively. First, we prove the inequalities for n = 0. 

/ rp / \l/(l-m) 

( t -At j M °' ( T 2 -At ) wm ^ e su P er solutions of (13. ip . (13.21) if 



m / T \ m /( 1 - m ) im / T \ 1 /( 1 ~ m ) / T \ 1 /( 1 ~ m ) 



1 - m V % - At 7 u 1 - m V To - At 7 u V To - At 



m f To \ m /(i-m) m f To \ 1/(1_m) / T 9 \ VC 1 -" 1 ) 



<+T— r <+At (iH-artio) > 



1 - m V To - At 7 u 1 - m V To - At 7 u V To - At 



that is, {1 _Z)t 2 u o 1 + -av >0 and (1 _^ )Ta ug* + Aw - aw > 0. 
From (15.21) . these inequalities are satisfied. 

Assume now that the estimates are satisfied at the order n.We prove it is satisfied at the order 

71+1. 

/ \ l/(l-m) / \ l/(l-m) 

The functions f^f^ 2 -) «„, ( jfz^J w « wi U be supersolutions of (iO) . (I3T2I) if 

771 I 777. I 

— u™ + Au n - av n > 0, — v™ + Av n - au n > 0. 

l — ml 2 — t n l — ml 2 — t n 

Since u n , v n are the solutions at the time level t n , we get : 

+ v4u n - aw n = — — — u n — ■ u n u 



™ ( i,m T 2 tn-1 _ m-1 

1 - mT 2 - t n (1 - m) At \ T 2 - t„ 
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rn 



1-mTo-L 



+ Av, n 



m 



,T 2 -t 



au r 



n-l 



m)At 



To - t r . 



By using the recurrence hypothesis, we obtain the estimates (15. 3p . 



Remark 5.7. // Au — av > and Av — au > 0, then T 2 = +oo and we obtain: u n+ \ < 

Urn ^n+l ^ V n . 



From these two lemmas, we obtain the inequalities 

~ < «r - W < 

and analogous inequalities for v n . Then, we obtain the convergence of the scheme as in [3]. 
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